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ABSTRACT 


An  analysis  of  the  heterogeneous  ignition  of  a 
solid,  under  convective  heating  by  a hot  gas,  for  large 
values  of  the  activation  energy  reveals  the  existence 
of  two  ignition  regimes.  In  a first  regime,  for  large 
values  of  the  frequency  factor,  a short  reactive  stage 
ending  in  thermal  runaway  follows  a long  inert  stage. 
For  small  values  of  the  frequency  factor  the  chemical 
heat  release  can  be  neglected  in  a first  stage  but  not 
in  a much  longer  second  stage  when  the  surface  temper- 
ature is  close  to  the  hot  gas  temperature;  a critical 
value  is  found  for  the  frequency  factor  below  which  no 
thermal  runaway  or  well  defined  ignition  time  exists. 


Closed  form  expressions  for  the  ignition  time 
are  given  for  both  regimes. 

D D C 
7?  JZtSjtETol 

Ef)  CCD  YC 


DISTRIBUTION  STATEMENT 

Approved  ten  public  release; 
Distribution  Uc.lim.itad 


FEB  16  19?? 

D 


I INTRODUCTION 


This  paper  deals  with  the  analytical  description 
of  the  process  of  heterogeneous  ignition  of  a solid 
under  convective  heating  by  a hot  gas.  The  heat  exchange 
between  the  gas  and  the  solid  is  modelled  by  a Newtonian 
law,  so  that  the  heat  flux  from  the  gas  to  the  solid 
is  written  as  proportional  to  the  difference  in  temper- 
ature of  the  gas  and  that  of  the  surface  of  the  solid. 
The  exothermic  heterogeneous  reaction  is  considered 
to  follow  Arrhenius  kinetics  with  an  activation  energy 
large  compared  with  the  thermal  energy  of  the  solid. 

The  effects  of  interphase  mass  transfer  are  neglected 
here . 


For  a review  of  previous  numerical  and  analytical 
work  related  to  this  model,  see  the  review  by  Merzhanov 
and  Averson^ . An  asymptotic  analysis  for  large  activation 
energies  of  heterogeneous  ignition  under  conductive 
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heating  by  a hot  gas  was  presented  by  Lin an  and  Crespo 
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and  by  Kindelan  and  Linan  , who  took  into  account  the 
effects  of  gasification  accompanying  the  heterogeneous 
reaction.  An  asymptotic  analysis  of  ignition  under  con- 
vective heating  by  a hot  gas,  when  an  exothermic  reaction 
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An  analysis  of  the  heterogenous  ignition  of  a solid,  under  convective  heating  by 
a hot  gas,  for  large  values  of  the  activation  energy  reveals  the  existence  of  twc 
ignition  regimes.  In  a first  regime,  for  large  values  of  the  frequencies  factor, 
i short  reactive  stage  ending  in  thermal  runaway  follows  a long  inert  stage.  For 
small  values  of  the  frequency  factor  the  chemical  heat  release  can  be  neglected 
in  a first  stage  but  not  in  a much  longer  second  stage  when  the  surface  temper  i- 
ture  is  close  to  the  hot  gas  temperature;  a critical  value  is  found  lor  t h» 
frequency  factor  below  which  no  thermal  runaway  or  well  defined  ignition  um, 
exists.  Closed  form  expressions  for  the  ignition  1 1 me  are  given  f.i  t>.  t h i •>•.*»  » 
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occurs  in  the  condensed  phase  has  been  recently  completed 
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by  Niiooka  and  Williams  . The  analysis  presented  here 
for  heterogeneous  ignition  is  parallel  to  that  of  Niiooka 
and  Williams. 

The  chemical  heat  release,  which  is  a very  sensi- 
tive function  of  the  surface  temperature,  becomes  com- 
parable to  the  convective  heating  only  when  the  surface 
temperature  approaches  a value  Tc>  In  Section  III  we 
consider  the  cases  when  the  frequency  factor  is  large 
enough  so  that  Tc  is  intermediate  between  the  initial 
solid  temperature  and  the  gas  temperature.  In  this  case 
there  is  a first  long  stage  of  inert  heating  followed 
by  a short  transition  stage,  when  the  chemical  and  con- 
vective heating  are  comparable,  ending  in  a precipituous 
rise  in  temperature  or  thermal  runaway  at  a well  defined 
ignition  time.  In  Section  IV  an  analysis  of  the  ignition 
or  heterogeneous  reaction  process  is  presented  for  smaller 
values  of  the  frequency  factor  when  Tc  becomes  close 
to  the  hot  gas  temperature.  A critical  value  is  found 
for  the  frequency  factor  below  which  no  thermal  runaway 


occurs  . 
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II  FORMULATION 


When  a semiinfinite  body  is  heated  at  its  surface 
by  a flux  proportional  to  the  difference  in  its  surface 
temperature  and  that  of  the  source,  the  equations  des- 
cribing the  temperature  distribution  can  be  written  in 
non-dimensional  form  as  the  heat  conduction  equation, 


9_e 

9 t 


9^0 


9x‘ 


T 0 


x>  0 , 


t>0 


(1) 


to  be  solved  with  the  initial  condition 


0 (x , 0) =0  ^ 


which  is  also  the  temperature  of  the  solid  gas  from 
the  surface  at  all  times.  In  addition  we  have  the  sur- 
face condition 


x = 0 , t.-0  : - |^-=l-0+Be"0a/G  (3) 

where  we  have  assumed  that  an  exothermic  heterogeneous 
reaction  with  Arrhenius  kinetics  takes  place  at  the 


surface  of  the  solid. 


The  temperature  has  been  made  non-dimensional 
with  the  temperature  of  the  heat  source;  G is  the 

Q. 

non-dimensional  activation  temperature.  The  character- 
istic length  and  time  used  in  the  nondimensional  length 
x,  time  t and  frequency  factor  B have  been  obtained 
from  the  thermal  diffusivity,  heat  conduction  coefficient 
and  overall  heat  transfer  coefficient. 

As  is  well  known  there  is  a relation 


0-0  . 
1 


r0  l-e+Be"0^0 

J° — 7FF — J 

v t-t. 


(4) 


between  the  surface  temperature  and  heat  transfer  flux 
at  the  surface. 

We  shall  give  in  this  paper  a description  of  the 
solution  of  Eq(4)  for  large  values  of  0a,  first,  for 
fast  reactions,  corresponding  to  large  values  of 
B exp(-Oa),  in  Section  III,  and,  second,  for  slow  reac- 
tions, corresponding  to  values  of  0aB  exp(-0a)  or  order 
unity  or  small  in  Section  IV. 


Ill  IGNITION  WITH  FAST  REACTIONS 


For  B=0,  or  for  values  of  B and  0 such  that 
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^ 0 / 0 

Be  a <<(1-0),  the  surface  temperature  history  is  given 
by  0=0^(t),  the  inert  temperature  value  given  by  the 


solution  of 


t i-OjCtj) 


(t)-e.-  -i-  /- 

1 /7  0 /f^f; 


(5) 


takes  the  value 


( 1 ~ 0 £ ) / ( 1 -0 = eterfc/F 


(6) 


For  large  values  of  the  nondimensional  activation  tem- 
perature 0 and  values  of  B such  that 
a 


l-0c  = Bexp(-0a/0c)  , ei<Qc<1  (7) 


There  is  an  initial  inert  stage,  as  outlined  by  Lin  an 
and  Williams  , where  the  temperature  0 takes  the 
inert  value  of  Eq(6),  except  for  exponentially  small 
terms  ( ~ exp ( 0 /0  -0  / 0 ) ; the  inert  stage  is  followed 

3 C a 

by  a short  transition  stage  around  t=t  , where  t is 

c c 

given  by 
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Tht  transition  stage  ends  with  a thermal  runaway  at 
a well  defined  ignition  time  t^,  given  gy  the  relation 


n -0.431, 
B = e (- 


8li 


1/2 


e e ' 

a f i 


ea/0f  i 


(9) 


where 


ofi.ef(t.) 


and 


d0 


’fi  = 


d t 


•(ti) 


as  it  can  be  shown  through  a procedure  paralleling  that 

2 

used  by  Lilian  and  Crespo  . 

Notice  that  the  characteristic  time  of  the  tran- 
sition stage  is  of  order  1/0  B'(t  ),  and  becomes  very 
large  when  0 -*1  because  0^(tc)->-O;  in  this  case  the  analy- 
sis leading  to  Eq(9)  for  the  ignition  time  fails.  We  shall 
here  present  an  analysis  of  the  solution  of  the  system  of 
Eqs(l)-(3)  for  values  of  B such  that  the  relation  (7)  is 
satisfied  only  for  values  of  0 >1  , In  fact  there  is  a 
distinguished  limiting  case  corresponding  to  values  of 
very  close  to  1 , which  we  shall  analyze  in  detail  below. 


IV  IGNITION  WITH  SLOW  REACTIONS 


We  shall  describe  here  the  solution  of  Eq(4)  for 


; 

large  values  of  and  values  of  B such  that 


6=0  Bexp(-0  ) 
S 3 


00) 


is  of  order  unity.  In  this  case  the  temperature  distri- 
bution will  be  given  by  Eq(6)  until  t is  large  enough 
so  that 


1-0,.=  (1  — 0.)// ii~t 
f i 


Small  temperature  increments,  or  order  1/0  , around 

a 

0=1  suffice  to  multiply  the  reaction  term  by  a factor 
e making  it  comparable  or  more  important  that  the  new- 
tonian  heating  term.  We  shall  write  Eq(4)  in  terms  of 
the  variable. 


©=*/0a+°f 


(11) 


and  notice  that  in  first  approximation  we  can  write 
the  Arrhenius  exponent  in  the  form 


-0a/0 

e 


- 0 
e 


ac^-0a(l-0i) /y'Vt 


(12) 


which  is  the  correct  asymptotic  form  of  the  exponent 
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when  the  reaction  term  can  not  be  neglected.  With  this 
approximation  Eq(4)  becomes 


♦ - 


~4i  + 6 e 


-a/  v'Ti+'l; 


Z7 


-dt, 


/ t-t 


where 


(13) 


a=  0 (l-0.)/ZiT 
a x 


04) 


The  integral  equation(l 3)  , which  is  of  the  Abel 
type,  will  be  solved  below  for  large  values  of  a in 
the  distinguished  limiting  case  when  6 is  of  order  1 . 
To  obtain  the  asymptotic  solution  we  begin  by  writing 
it  in  terms  of  the  variable 


o=t/ a2 


(15) 


to  obtain 


i 0 ,/  uj  if)  (oi ) ”1  / 

_J_  r _uKqi)+*e  —ioi 

0 z^-oT 


which  can  be  inverted  to  yield 


i)j  - 1 / Z^_  1 


d^  ( 0 1 ) 


/ ■ 

ct  Z"  7t  Za-ai 


(16) 
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We  then  introduce  the  expansion 


\p  = \L  + — ^ — ill  + 

O O 1 


(18) 


in  Eq(17)  to  obtain 


i - 1 / /cT 

\b  e T o = o e 
T o 


(19) 


and 


a d\ |»  (o) 

f (♦  -I)-  -i-  f 2 

> 0 ,V  ° 


(20) 


Eq  (1 9)  t giving  the  first  approximation  for  \p(a), 
has  two  branches;  the  one  with  physical  meaning  is  the 
lower  one  and  it  shows  i|>  growing  with  afrom  0 to  an 
asymptotic  value  given  by 


\fiaexp(-\)»a)=6 


(21) 


if  <5  is  smaller  than  a critical  value,  6 c = l / e . For  super- 


This  equation  could  also  be  obtained  from  Eq(3),  written 
in  terms  of  the  variable  f , by  neglecting  the  heat 
conducted  df/dx  to  the  solid  and  by  linearizing  the 
Arrhenius  exponent. 





Udi  M.  > I— i.MOM  m 


PNHP 
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supercritical  values  of  6,  <5>  1 / e , ip  grows  from  0 to 


1 when  a grows  from  0 to  a critical  value  given  by 


6exp(-l/ /o  ) = 1 / e 


(22) 


In  this  supercritical  solution  when  a approaches  o^, 

t|j  - t because  t/; 1 and  d^Q/do-»-«>,  The  expansion  in  Eq(18) 

fails  and  should  be  replaced  by  another  involving  the 

3/2 

stretched  variable  (o-oc)w  • The  first  approximation 

for  i \>  in  this  transition  stage  shows  a thermal  runaway 

at  an  ignition  time  diferent  from  a ^ by  an  amount  or 
-3/2 

order  a . Thus  Eq(22)  can  be  used  in  first  approxi- 
mation to  calculate  the  ignition  time  for  supercritical 
values  of  6.  In  terms  of  the  original  variables  the 
ignition  time  t^  is  given  by  the  relation 


0 (l-0.)//ut.  = ln{e0  Be  a) 
ail  a 


(23) 


For  subcritical  values  of  <$ , that  is  for  values  of  B 


smaller  than  a critical  Bc> 


B e = l/e0 
c a 


(24) 


no  ignition  time,  leading  to  high  temperatures  exists, 


-li- 


the chemical  reaction  proceeds  at  an  asymptotic  temper- 
ature given  by  Eq(21)  which  is  close  to  the  source  tem- 
perature. 

V CONCLUSIONS 

A closed  form  relation  has  been  found  between  the 
nondimensional  frequency  factor  B and  the  ignition  time 
t^  by  solving,  for  large  values  of  the  nondimensional 
activation  energy,'  the  equations  (l)-(3),  or  the  integral 
equation(4),  which  describe  the  heterogeneous  ignition 
under  convective  heating  by  a hot  gas. 

The  analysis  sketched  in  Section  III  led  to  the 
closed  form  relation(9)  for  the  ignition  time  for  values 
of  B such  that  relation(7)  is  satisfied.  The  analysis 
of  Section  IV  corresponds  to  values  of  B such  that  6, 
given  by  Eq(10),  is  of  order  unity;  the  analysis  indi- 
cated the  existence  of  a short  inert  stage  where  the 
surface  temperature  is  given  by  Eq(6),  followed  by  a 
long  reactive  stage  when  the  surface  temperature  is 
given  in  first  approximation  by 

0 = 1 + 0 (i|j  -l//o) 

a 0 


(25) 


-12- 


and  the  nondimensional  increment  in  surface  temperature 
ty  due  to  the  reaction,  is  given  by  Kq(19),  with 

c = TTt/02  (1-0  .)  2 (26) 

a 1 

For  values  of  6 below  a critical  value  6c  = l/e,  ^ grows 
asymptotically  to  a value  i|i  given  by  Eq(21),  while 

ci 

for  larger  values  of  6 a thermal  runaway  is  found  at 
in  ignition  time  t^  given  by  Eq(23). 
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